PROBING BEYOND THE EPOCH OF HYDROGEN REIONIZATION WITH 21 CENTIMETER 
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ABSTRACT 

We use numerical simulations of hydrogen reionization by stellar sources in the context of ACDM 
cosmogonies to investigate the 21 (1 + z) cm radio signal expected from the diffuse intergalactic medium 
(IGM) prior to the epoch of reionization breakthrough at redshift z- lon . Two reionization scenarios are 
analyzed in details: an 'early reionization' case with z- lcm « 13, consistent with the recent discovery by 
the WMAP satellite of a large optical depth to Thomson scattering, and a 'late reionization' case with 
zion ~ 8. It is a generic prediction of these models that the background of Lya photons produced by 
the early generation of stars which ultimately ionize the universe will be sufficiently intense to make 
intergalactic neutral hydrogen visible against the cosmic microwave background during the 'grey age', 
i.e. Zion z < 20. Depending on the redshift of reionization breakthrough, broad-beam observations at 
frequencies < 150 MHz (below 100 MHz for the 'early reionization' scenario) with the next generation of 
radio telescopes should reveal angular fluctuations in the sky brightness temperature in the range 5-20 
mK (lcr) on scales below 5 arcmin. 

Subject headings: cosmology: theory - diffuse radiation - galaxies: evolution - intergalactic medium 
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1. INTRODUCTION 

Despite much recent theoretical and observational 
progress in our understanding of the formation of early 
cosmic structures and the high redshift universe, many 
fundamental questions remain, at best, only partially an- 
swered. We know that at least some galaxies and quasars 
were already shining when the universe was less than 10 9 
yr old. But when did the first luminous objects form, 
and what was their impact on the surrounding intergalac- 
tic gas? While the excess H i absorption measured in the 
spectra of z ~ 6 quasars in the Sloan Digital Sky Sur- 
vey has been interpreted as the signature of the trailing 
edge of the cosmic reionization epoch (Becker et al. 2001; 
Fan et al. 2002), the recent analysis of the first year data 
from the Wilkinson Microwave Anisotropy Probe ( WMAP) 
satellite infers a mean optical depth to Thomson scatter- 
ing r e ~ 0.17, suggesting that the universe was reionizcd 
at higher redshift (Kogut et al. 2003; Spergel et al. 2003). 
This is an indication of significant star-formation activ- 
ity at very early times. In currently popular cosmological 
scenarios it is an early generation of subgalactic stellar 
systems - aided perhaps by a seed population of accret- 
ing black holes in their nuclei - that may have generated 
the ultraviolet radiation that ended the cosmic 'dark ages' 
and reionized most of the hydrogen in the universe around 
redshift 8-15. The detailed thermal history of the baryons 
during and soon after these crucial formative stages de- 
pends on the power-spectrum of density fluctuations on 
small scales and on a complex network of poorly under- 
stood 'feedback' mechanisms, and remains one of the miss- 
ing link in galaxy formation and evolution studies. 

An alternative way to probe the end of the cosmic 'dark 
ages' and discriminate between different reionization his- 
tories is through 21 cm tomography. It has also long been 
known (e.g. Field 1959; Sunyaev & Zel'dovich 1975; Hogan 
& Rees 1979; Scott & Rees 1990; Subramanian & Pad- 
manabhan 1993) that neutral hydrogen in the intergalactic 



medium (IGM) and gravitationally collapsed systems may 
be directly detectable in emission or absorption against 
the cosmic microwave background radiation (CMB) at the 
frequency corresponding to the redshifted 21 cm line (asso- 
ciated with the spin-flip transition from the triplet to the 
singlet state). Madau, Meiksin, & Rees (1997, hereafter 
MMR) first showed that 21 cm tomography could pro- 
vide a direct probe of the era of cosmological reionization 
and reheating. In general, 21 cm spectral features will dis- 
play angular structure as well as structure in redshift space 
due to inhomogeneities in the gas density field, hydrogen 
ionized fraction, and spin temperature. Several different 
signatures have been investigated in the recent literature: 
(a) the fluctuations in the 21 cm line emission induced by 
the 'cosmic web' that develops at early times in cold dark 
matter (CDM)-dominated cosmologies (Tozzi et al. 2000) 
and by 'minihalos' with virial temperatures below 10 4 K 
(Iliev et al. 2002,2003); (b) the global feature ('reioniza- 
tion step') in the continuum spectrum of the radio sky 
that may mark the abrupt overlapping phase of individual 
intergalactic H n regions (Shaver et al. 1999); (c) and the 
21 cm narrow lines generated in absorption against very 
high redshift radio sources by the neutral IGM (Carilli, 
Gnedin, & Owen 2002) and by intervening minihalos and 
protogalactic disks (Furlanetto & Loeb 2002). 

In this paper, we use numerical simulations of hydrogen 
reionization by stellar sources in the context of ACDM 
cosmogonies to estimate the radio signal expected from 
the diffuse, low density IGM during the 'gray age' (Car- 
illi et al. 2002). This is the era between the formation 
of the first luminous objects and the epoch of complete 
reionization, when the IGM is a mixture of neutral, par- 
tially ionized, and fully ionized structures. The important 
point is that the early generation of stars that ultimately 
ionizes the universe will also generate a background radia- 
tion field of Lya photons arising from their collective red- 
shifted UV continua. It is the presence of a sufficient flux of 
Lya photons that renders the neutral IGM 'visible' against 
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the CMB by mixing the hyperfine levels ('Lya pumping'). 
Radio maps will show a patchwork (both in angle and in 
frequency) of emission signals from H I zones modulated 
by H ii regions where no signal is detectable against the 
CMB. 

The remainder of this paper is organized as follows. In 
§2 we summarize the physics of 21 cm excitation. In §3 
we show that, since a 'typical' early stellar population will 
radiate a large number of UV continuum photons with 
energies between 10.2 and 13.6 eV (hereafter Lya-c ) to 
which the IGM is transparent, there will be an extended 
period of time prior to full reionization where fluctuations 
in the gas density and hydrogen ionized fraction may be 
detected at 21 (1 + z) cm. Numerical simulations of hy- 
drogen reionization by Population III stars are presented 
in §4. Simulated radio maps are generated and the de- 
tectability of fluctuations in the sky brightness discussed 
in § 5. Finally, in § 6 we give our conclusions. 

2. LYMAN-ALPHA PUMPING 

A quick summary of the physics of 21 cm radiation will 
illustrate the basic ideas behind this work (see MMR for 
more details). The emission or absorption of 21 cm pho- 
tons from a neutral IGM is governed by the hydrogen spin 
temperature, Tg, defined by 



ni/no = 3exp(-T st /T s ), 



(1) 



where no and n\ are the number densities of atoms in 
the singlet and triplet n = 1 hyperfine levels and T* = 
0.068 K is the excitation temperature of the transition. 
In the presence of only the CMB radiation with Tomb = 
2.73 (1 + z) K, the spin states will reach thermal equilib- 
rium with the CMB on a timescale of T*/(Tcmb^4io) ~ 
3 x 10 5 (1 + z) -1 yr (A w = 2.85 x lO^s" 1 is the spon- 
taneous decay rate of the hyperfine transition of atomic 
hydrogen), and intergalactic H i will produce neither an 
absorption nor an emission signature. Thus, a mechanism 
is required that decouples Tg and Tcmb, e.g. by coupling 
the spin temperature instead to the kinetic temperature 
Tk of the gas itself. Two mechanisms are available, spin- 
exchange collisions between hydrogen atoms (Purcell & 
Field 1956) and scattering by Lya photons (Wouthuy- 
sen 1952; Field 1958; Field 1959). The first process 
proceeds at a rate that is too small for realistic IGM 
densities at the rcdshifts of interest, although collisions 
will be important in dense regions with gas overdensities 
5 > 20[(1 + z)/10]~ 2 , like virialized minihalos (Iliev et al. 
2002,2003; Furlanetto & Loeb 2002). It is Lya pumping 
that dominates instead in the diffuse IGM, by mixing the 
hyperfine levels of neutral hydrogen in its ground state via 
intermediate transitions to the 2p state, the Wouthuysen- 
Field process. In this case, 



P n = 



47rJ a cr 
hpv a 



(2) 



is the rate at which Lya photons of energy hpu a in an 
isotropic background of intensity J a are scattered by an H 
atom in the gas, and 



Pth 



27AioTcmb 
471 



7.6 x 10~ 13 (1 + z) s" 



(3) 



is the critical ('thermalization') scattering rate which, if 
much exceeded, would drive Tg away from Tcmb and to- 
wards Tk (MMR). In the above equations a a = J a v dv = 
Tre 2 f a /m e c = 0.011 cm 2 Hz is the integrated Lya scatter- 
ing cross section. The condition P a > P t h can also be 
rewritten as 

J a > Jth = 9xl0~ 23 (l + z) ergscm^s^Hz^sr- 1 . (4) 

At z = 8 this thermalization rate corresponds to 
Pt\i v al \cn^a a ) < 1/2 Lya photons per hydrogen atom, 
where 7jh is the intergalactic mean hydrogen density. 

3. UV SPECTRA OF EARLY STELLAR POPULATIONS 

At the epochs of interest, Lya pumping will efficiently 
decouple T5 from Tqmb if J a ^ 10 -21 ergs cm~ 2 s _1 Hz -1 
sr _1 . The issue is then whether this critical value of J a 
will be reached well before the epoch of complete reion- 
ization, when large regions of the IGM are still neutral. 
The process of reionization begins as individual sources 
drive isolated H 11 regions in the surrounding IGM and 
ends when these ionized bubbles overlap and fill the inter- 
galactic volume. The volume filling factor of H 11 regions 
in a clumpy medium reaches unity when about 1 + t/t Tec 
ionizing photons have been emitted per hydrogen atom, 
where i rcc is the volume-averaged hydrogen recombina- 
tion timescale (Madau, Haardt, & Rees 1999). Based on 
the density field in large scale cosmological simulations, 
only a few Lyman-continuum (Lyc) photons per hydro- 
gen atom may be needed (Gnedin 2000; Miralda-Escude, 
Hachnelt, & Rees 2000) to keep the gas in mildly overdense 
regions and intergalactic filaments ionized against radia- 
tive recombinations at z ~ 7; many photons, however, 
will not escape local absorption, i.e. only a small frac- 
tion, ,/csc, of the Lyc radiation emitted by massive stars is 
expected to escape the dense sites of star formation into 
galaxy halos and the intergalactic space. Moreover, at 
the earliest epochs of structure formation in CDM cos- 
mologies the smallest nonlinear objects are the numerous 
halos that condense with masses just above the cosmo- 
logical Jeans mass. Such 'minihalos' are not yet fully re- 
solved nor is the process of their photoevaporation cap- 
tured in large-scale three-dimensional cosmological simu- 
lations. This sub-grid structure makes the treatment of 
IGM clumping in all presently available simulations only 
approximate (Haiman, Abel, & Madau 2001). The error 
in small-scale clumping (that leads to underestimating the 
total recombination rate) can be partly compensated by 
the assumption of a small escape fraction. 

The early generation of massive stars likely responsi- 
ble for reionization will also generate a background ra- 
diation field of photons with energies between 10.2-13.6 
cV to which the IGM is transparent. As each of these 
Lya-c photons gets redshifted, it will ultimately reach the 
Lya transition energy of 10.2 eV, scatter resonantly off 
neutral hydrogen, and mix the hyperfine levels. The ob- 
servability of the pre-reionization IGM depends therefore 
on the UV spectrum of the firsts stars, i.e. on the num- 
ber of Lya-c photons emitted per H-ionizing photon, and 
on the escape fraction of Lyc photons into the IGM. Fig- 
ure 1 shows the number of continuum photons per octave 
near the Lya frequency emitted by a stellar population, 
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N a . The UV spectral energy distribution is character- 
ized by a strong Lyc break the size of which depends on 
age, initial mass function (IMF), metallicity, and star for- 
mation history. We have therefore normalized N a to the 
total number of hydrogen-ionizing photons emitted, 7Vj on . 
The values shown have been computed assuming a sim- 
ple stellar population of metal-free 'Population III' stars: 
the simulations described in the next section use synthetic, 
time-dependent spectra of zero-metallicity stars to model 
the emission properties of stellar sources. For compari- 
son, the figure also depicts the number of Lya photons 
emitted by low metallicity stars (Z = Z & /50). Note how 
the hotter temperatures of Population III stars produce a 
smaller Lyc break hence a smaller N a /N[ on ratio. Also, 
after integrating over stellar age, a starburst is character- 
ized by a ratio N a /Ni on that is a few times higher than 
in the constant star formation rate case (this is because 
of the shorter lifetime of the massive stars that produce 
radiation above 1 ryd). It is possible that the starburst 
mode may be more relevant for the galaxies responsible 
for reionization, as stellar winds and supernovae can eas- 
ily expel the gas out of the shallow potential wells of low 
mass halos, after the first burst of star formation occurs 
(Madau, Ferrara, & Rees 2001). 
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Fig. 1. — Number of continuum photons per octave in frequency 
just shortward of Lya, N a , emitted as a function of age by a stel- 
lar population with a Salpeter IMF (dot-dashed lines) and a Scalo 
(long-dashed lines) IMF. Upper set of curves assumes an instanta- 
neous burst of star formation, lower set a constant star formation 
rate. We have normalized N a to the number of hydrogen-ionizing 
photons per octave emitted just shortward of the Lyman limit, N- lon . 
The population synthesis values have been computed for low metal- 
licities (Z = 0.02 Zq), and are based on an update of Bruzual & 
Chariot's (1993) libraries. Solid line: same for zero-metallicity Pop- 
ulation III stars with a Salpeter IMF. Short-dashed line: same for 
zero-metallicity Population III stars with a Larson IMF (see Ciardi 
et al. 2001 and references therein). 

As mentioned above, only an amount iVphot — /esc Mon 
of Lyc photons will actually be able to escape from the 
galaxies and break into the IGM. A smaller escape fraction 
will generate a stronger Lya diffuse radiation background 
for a given ionizing metagalactic flux. Taking a fiducial 
ratio N a /Ni on > 4 and a few H- ionizing photons per atom 
to keep the IGM fully photoionized, N p h t/Nn w 5 (see 
Fig. 7 in Ciardi, Stoehr, & White 2003, hereafter CSW), 



these simple estimates imply 

^>4%^/ o - c 1 > 200(/ csc /0.1)- 1 (5) 

Lya 'continuum' photons emitted per hydrogen atom close 
to the reionization epoch. This is more than two orders 
of magnitude larger than the thermalization rate for effi- 
cient level mixing by Lya pumping. It appears then to be 
a generic prediction of cosmological reionization scenarios 
by stellar sources that the Lya background from the early 
generation of stars which ultimately ionize the universe 
should be sufficiently intense to make intergalactic diffuse 
neutral hydrogen visible against the CMB well before the 
epoch of full reionization. As shown below, detailed sim- 
ulations of IGM reionization that include cosmological ra- 
diative transfer effects will confirm this basic result. 

4. NUMERICAL SIMULATIONS OF IGM REIONIZATION 

In this and the following section we use the numerical 
simulations of hydrogen reionization by stellar sources de- 
scribed in CSW and Ciardi, Ferrara, & White (2003, here- 
after CFW) to investigate the radio signal expected from 
the diffuse IGM prior to reionization. We give here a brief 
summary of the main features of the simulations, relevant 
to the present study. In the above papers, the reionization 
process is studied through a combination of high resolution 
N-body simulations (to describe the distribution of dark 
matter and diffuse gas, Springel, Yoshida, & White 2001; 
Yoshida, Sheth, & Diaferio 2001), a semi-analytic model 
of galaxy formation (to track gas cooling, star formation, 
and feedback from supernovae, Kauffmann et al. 1999; 
Springel ct al. 2001) and the Monte Carlo radiative trans- 
fer code CRASH (to follow the propagation of ionizing pho- 
tons into the IGM, Ciardi et al. 2001; Maselli, Ferrara, & 
Ciardi 2003). The simulations, spanning the redshift range 
z = 6 — 20, are performed for a ACDM cosmology with 
(Q m ,Q A) (T 8 , h, Ofe, n) = (0.3,0.7,0.9,0.7,0.04,1). The 
simulation box has a comoving length of L = 20/i _1 Mpc 
and a particle mass of M p = 1.7 x 10 8 /i _1 M Q . This choice 
of the parameters allows to simulate a region of the uni- 
verse with 'mean' properties, thus avoiding biases due to 
cosmic variance on small scales and, at the same time, to 
resolve the objects that produce the bulk of the ionizing 
radiation (Ciardi 2002). Each output of the simulation 
provides the gas distribution inside the box together with 
a galaxy catalog containing, among other quantities, their 
positions, stellar masses, and star formation rates. As dis- 
cussed in § 5, two sets of simulations will be used, with 
different choices for the galaxy emission properties. At 
each output, the propagation into the IGM of the ionizing 
photons emitted by the simulated galaxy population is fol- 
lowed with the radiative transfer code CRASH. The code has 
been tested to ensure a good convergence in the ionization 
fraction for a wide range of densities. The calculation is 
performed over aJV= 128 3 grid. The gas (pure hydrogen) 
distribution is assumed to track the dark matter. While 
radiation above 13.6 eV from each source is fully absorbed 
at the ionization front surrounding it, the Lya-c radia- 
tion is free to propagate and fill the space in between the 
H II regions. Since numerous sources are contributing to 
the Lya background throughout the universe, its flux is 
nearly isotropic in the cosmic frame of reference, and is 
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given by 



e(v',z') 



{l + zf 



dt 



dz 1 



dz , (6) 



(1 + z') 3 

where the proper emissivity e(v' , z') can be written as 
(1 + z'f 



L 3 



-M*{z')S{v'). 



(J) 



Here, v' = v a (l + z')/(\ + z), M^(z') is the total stellar 
mass in the simulation box at redshift z', and S(v') is the 
time-dependent spectral energy distribution at frequency 
v' . Since the UV photon production is rapidly decreasing 
with stellar age, M^(z') is the mass of newly (i.e. in the 
timescale between two outputs of the simulation) formed 
stars, i.e. we neglect the contribution from stars formed 
in previous outputs. In equation (6) we have assumed 
that Lya photons are not attenuated by the intervening 
IGM, and that the IGM is opaque to ionizing photons 
prior to reionization, z > z- lon . Thus, the integration over 
redshift is carried out only up to 1 + z max = 4(1 + z)/3. 
In Figure 2 the resulting Lya background intensity is de- 
picted as a function of redshift, together with the volume- 
averaged ionization fraction x v , for the two simulation 
runs. The figure clearly shows that the expected diffuse 
flux of Lya photons should be sufficiently intense to decou- 
ple the spin temperature from the CMB during the 'grey 
age' from redshift ~ 20 to complete reionization. 
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FlG. 2. — Solid lines: redshift evolution of the 'continuum' 
Lya background, J a , predicted by our simulations of cosmological 
reionization by stellar sources. Dashed lines: volume-averaged ion- 
ization fraction, x v . For each set of curves, the thick (thin) line refers 
to the 'early reionization' ('late reionization') scenario. Dotted line: 
Thermalization rate, J t ^, needed to efficiently decouple the spin tem- 
perature of neutral hydrogen from the CMB via Lya pumping (see 
text for details). 



5. RADIO EMISSION FROM THE DIFFUSE IGM 

A beam of 21 cm radiation passing through a neutral 
hydrogen patch causes absorption and induces emission. 
In the radio range the exponential in equation (1) is close 



to unity and stimulated emission cancels most of the ab- 
sorption. The 21cm absorption coefficient depends then on 
the spin temperature. If the patch is located at redshift 
Z, has (mean) spin temperature Ts, cosmic overdensity <5, 
neutral hydrogen fraction (1 — x), angular diameter on the 
sky larger than the beamwidth of the radio telescope, and 
a radial velocity width broader than the bandwidth, its 
optical depth at 21(1 + z) cm is given by (MMR) 



0.006 r 



T c 



MB 



(8) 



where we have included in the function T all dependences 
on the density and ionization state of the patch, and on 
cosmological parameters, 



0.02 



1 



10 



0.3 

S2m 



1/2 



(9) 



Here the Hubble constant is Ho = 100 h km s 1 Mpc 1 
and we have used the approximation H{z) rs \/fl m Ho(l + 
zf' 2 . 

Because of the weakness of the magnetic hyperfine tran- 
sition, r will typically be much less than unity. In the 
rest frame of the patch, the radiative transfer equation 
yields for the brightness temperature through the translu- 
cent IGM T b = T CM Be" T + T s (l - e~ r ). The flux ob- 
served at Earth can be expressed by the differential an- 
tenna temperature between this patch and the CMB, 
ST b = {T b -T CM B){l + z)- 1 or 



rrr, ~ ^CMB 

6T b ~ — r 

1 + 2 



0.016 K T 1 



TcMB -). (10) 



If P a > Pth and the IGM has been preheated by pri- 
mordial sources of radiation, T$ will be much larger than 
Tomb, an d the universe will be observable in 21 cm emis- 
sion at a level that is independent of the exact value of 
Ts- If preheating is negligible, the adiabatic expansion of 
the universe will lower the kinetic (hence spin) tempera- 
ture of the gas well below that of the CMB, and the IGM 
will be detectable in absorption against the CMB. The 
energetic demand for heating the IGM above the CMB 
temperature is meager, only ~ 0.004 eV per particle at 
z ~ 10. Consequently, even relatively inefficient heating 
mechanisms may be important warming sources as struc- 
ture develops before the universe was actually reionized, 
such as Lya heating (MMR) and X-ray heating (MMR; 
Vcnkatesan, Giroux, & Shull 2001) in low density regions, 
and shock heating of the gas in high density regions (Carilli 
et al. 2002). In the cosmological simulations of Carilli et al. 
(2002), the IGM is reionized 'late' but the volume-averaged 
kinetic and spin temperatures of the gas exceed signifi- 
cantly Tomb already at redshifts 14 and 10, respectively. 
In the following we will focus on the radio emission signal 
expected from our simulations in the limit Tg 3> Tcmb- 

5.1. 'Late reionization' model 

The first simulation describes a 'late reionization' sce- 
nario in which the galaxy emission properties are com- 
puted assuming a Salpeter IMF, a time-dependent spec- 
trum typical of metal-free stars, and an escape fraction 
/ osc = 5% (CSW; run S5 in CFW). The choice of these 





Fig. 3. — Slices through the simulation box in the 'late reionization' model. The nine panels show the differential antenna temperature 
log (5T b /K) at redshifts (from top to bottom and left to right) z = 13.5, 12.8, 12.0, 11.3, 10.6, 9.9, 9.3, 8.7, and 8.1. The maps refer to the 
emission produced by a slice of gas of comoving thickness L/N 1 / 3 = 220 kpc. 



parameters leads to a reionization epoch z- lon sa 8 (defined 
here as the redshift where the volume-averaged ionization 
fraction, x v , reaches unity). We use the expressions de- 
rived in equations (8)-(10) together with our simulations to 
account for the spatial distribution of structures and their 
different ionization levels during the reionization process. 
Maps of the antenna temperature can be thus constructed 
at each redshift based on equation (10), with T in equation 
(9) derived from the simulations along each line-of-sight. 
In Figure 3, the resulting 5Tb maps are shown for differ- 
ent redshifts. The maps refer to the emission produced 
by a slice of the simulation box of comoving thickness 
L/N 1 / 3 = 220 kpc. Highly ionized regions (the dark areas 
in the map) initially occupy only a small fraction of the 
volume. At redshift 12, when about 90% of the IGM is still 
neutral, 6Tb reaches values ~ 0.1 K in mildly overdense re- 
gions with 6 greater than a few: the antenna temperature 
drops rapidly with cosmic time as reionization proceeds. 
By redshift of 11 or so, several bubbles have already over- 
lapped as the first sources of ionizing photons are highly 
clustered. Finally, by z = 8.7 most of the volume has been 
photoionized with only a few dense H i patches surviving. 

The distribution of 5Tb values is shown in Figure 4. At 
z = 18.5, most of the IGM is still neutral, and the bright- 
ness temperatures reflect density inhomogeneities still in 
the linear regime. As reionization gets under way, the 
distribution becomes double peaked due to the two-phase 
structure of the IGM: the rightmost peak corresponds 
to the neutral phase, while the left peak is associated 
with ionized gas. The pixel fraction corresponding to 
6Tb <C 10 -8 K is not drawn in the figure: this increases 
by two orders of magnitude as z- lon is approached. With 
decreasing redshift, more pixels are transferred from the 



right to the left peak. While the average 21 cm signal is 
two orders of magnitude lower than that of the CMB (and 
will be swamped by the much stronger non-thermal back- 
grounds that dominate the radio sky at meter wavelengths 
and that must be removed) , its fluctuations - induced now 
both by inhomogeneities in the gas density and in the hy- 
drogen ionized fraction - will greatly exceed those of the 
CMB. Brightness temperature fluctuations will be present 
both in frequency and in angle across the sky, and should 
be much easier to detect. 
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Fig. 4. — Distribution in our 'late reionization' simulation box of 
the differential antenna temperature 5T b (K) at redshifts z = 18.5 
(solid line), 12.8 (dotted line), 9.9 (small- dashed line), 8.7 (long- 
dashed line) and 8.1 (dotted- dashed line). 
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Fig. 5. — Upper panel: Expected rms brightness temperature 
fluctuations, (5T^ ) 1 / 2 , in our 'late reionization' simulation box as a 
function of beam size AO. A fixed bandwidth Av = 1 MHz has been 
assumed. Every curve corresponds to a different emission redshift or, 
equivalently, observed frequency f (, s . Lower panel: same in a test 
run with radiation sources switched off. Note that the box dimension 
is about 10 arcmin. 



Figure 5 (upper panel) shows the expected rms tempera- 
ture fluctuations relative to the mean, (ST 2 ) 1 / 2 , produced 
by a gas volume corresponding to a given bandwidth, 
Av = v b s Az/ (1 + z), and angular size A9, as a function of 
AO. The bandwidth is fixed for each redshift and equal to 
Av = 1 MHz. The signal peaks at z ~ 10.5 (~ 120 MHz), 
corresponding to the epoch when several high density neu- 
tral regions are still present, but H n occupies roughly half 
of the volume. As the mass variance is larger on small 
scales, at a fixed bandwidth (5T 2 ) 1 / 2 increases with de- 
creasing angular scale from about 5 to 20 mK. 
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Fig. 6. — Same as the upper panel of Fig. 5 for two different 
bandwidths, At/ = 0.1 (dotted lines) and 1 MHz (solid). 



At lower redshift (5T 2 ) 1 / 2 decreases and it drops dra- 
matically at z < 8.5, when the IGM gets close to complete 
reionization. Above redshift 11, the rms decreases again 
as the development of structures in the universe is less ad- 
vanced. Although inhomogeneities in the ionized fraction 
drop to zero as the IGM is mostly neutral, a patchwork 
of radio emission is produced by regions with overdensi- 
ties still in the linear regime and by voids. Note that the 
errors due to the box finiteness allow a correct evaluation 
of (5T 2 ) 1 / 2 only for AO < 8 arcmin. To better gauge the 
suppression effect of cosmological reionization on the sky 
brightness at 21 (1 + z) cm, we have plotted in the lower 
panel the rms temperature fluctuations expected in the 
same simulation but with radiation sources switched off. 
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Fig. 7. — Expected rms brightness temperature fluctuations in the 
'late reionization' case as a function of the observed frequency, v obs , 
for different choices of the beam angle, Ad. Notation is the same as 
in Fig. 5. The bandwidth is Av = 1 MHz. 



We have also checked that a smaller bandwidth of 
Av — 0.1 MHz gives roughly the same rms at high red- 
shifts, while this becomes larger on small scales as z de- 
creases (Fig. 6). This is related to the fact that inhomo- 
geneities in the gas density and in the hydrogen ionized 
fractions, which increase with decreasing redshift, are bet- 
ter resolved by a smaller bandwidth. The expected rms 
brightness temperature fluctuations as a function of the 
observed frequency, v bs, are shown in Figure 7 for our 
'late reionization' run, for different choices of beam size, 
and for a bandwidth Av = 1 MHz. The temperature fluc- 
tuations remain approximately constant over a wide range 
of redshift, drastically declining only when the IGM be- 
comes highly ionized. 

5.2. 'Early reionization' model 

Recently, the discovery by the WMAP satellite of a large 
optical depth to Thomson scattering (Kogut et al. 2003; 
Spergel et al. 2003) suggests that the universe was reion- 
ized at Zi on > 10. We have thus repeated the same analyses 
of § 5.1 for a simulation in which the galaxy emission prop- 
erties are computed assuming instead a Larson IMF (i.e. a 
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Salpeter function with a cutoff below 5 M©) and an escape 
fraction / osc = 20% (run L20 in CFW). These parame- 
ters lead to an 'early reionization' epoch, z- mn w 1 3, more 
consistent with the WMAP observations (see CFW for an 
extensive discussion). The overall, qualitative results are 
the same described in the previous Section for the 'late 
reionization' case. Here, we limit our discussion to the 
expected value of the rms brightness temperature fluctua- 
tions, shown in Figure 8. The same curves in Figure 8 and 
Figure 5 (upper panel) correspond to redshifts of equiva- 
lent volume averaged ionization fraction, e.g., x v ~ 0.5 is 
reached at z = 14.9 (10.4) in the 'early (late) reionization' 
case and it corresponds to the short-dashed line in both 
Figures. The behavior is similar to the 'late reionization' 
model, except that now the neutral IGM is detectable only 
at radio frequencies below 100 MHz. 



10 



A 













&v=l MHz 






. . . E 


z 


"„b. \ 




- 12.8 


103 MHz 




. 13.0 


100 MHz 




- 13.7 


97 MHz 




; 14.9 


90 MHz 


N = 


- 18.5 


73 MHz 


s - 

S - 



0.1 1 



A0(arcmin) 



Fig. 8. — Same as in Fig. 5 (upper panel) but in the 'early reion- 
ization' case. Note how rms brightness temperature fluctuations on 
all scales now drop below 2 mK at z < 13, making the IGM effec- 
tively undetectable in 21 cm radiation above radio frequencies of 100 
MHz. 



6. DISCUSSION AND CONCLUSIONS 

In this paper we have used numerical simulations of hy- 
drogen reionization by stellar sources in the context of a 
ACDM cosmogony to investigate the 21 (1 + z) cm radio 
background expected prior to the epoch of cosmological 
reionization. Here, in addition to the emission induced 
by the 'cosmic web' (Tozzi et al. 2000), we can also esti- 
mate the contribution from inhomogeneities in the hydro- 
gen ionized fraction. A further contribution comes from 
the minihalos with virial temperatures below the threshold 
for atomic cooling that form in abundance at high redshift, 
and are sufficiently hot and dense to emit collisionally- 
excited 21 cm radiation (Iliev et al. 2002,2003). These 
halos are not resolved in our numerical simulations. While 
the combined signal from many minihalos within a suffi- 
ciently large volume is in principle detectable, the 21 (1+z) 
cm flux from these objects will typically be dominated 
by the radio background from the neutral IGM (since the 
H I mass will always be larger in the non-collapsed phase; 
Oh & Mack 2003), except in the very earliest stages of 
structure formation when J a < J t h, i.e. at z > 20. At 



z > 12, the predicted brightness temperature fluctuations 
from minihalos on arcmin scales are typically much smaller 
than those expected from the diffuse IGM. As reionization 
proceeds, feedback effects (like photoevaporation) on mini- 
halos may also decrease their expected 21 cm signal. Here, 
we have focused on emission from the diffuse, low density 
gas. The search at 21 cm for the epoch of first light has be- 
come one of the main science drivers of the LOw Frequency 
ARray (LOFAR). Radio 21 cm tomography may probe the 
topology of reionization, map the cosmic 'gray age', i.e. 
the transition from the post-recombination universe to one 
populated with radiation sources, and effectively open up 
much of the universe to a direct study of the reheating and 
reionization epochs. On scales ~ 5 arcmin, predicted rms 
brightness temperature fluctuations can exceed 5 mK for 
a 1 MHz bandwidth. The expected sensitivity of LOFAR 
on these scales is about 10 mK with a 1000 h integration 
(for a confidence level of 5 times the noise and Av = 1 
MHz, see http://www.astron.nl/lofar/science/), enough to 
detect all fluctuations above 2 a. While remaining an ex- 
tremely challenging project due to foreground contamina- 
tion from unresolved extragalactic radio sources (Di Mat- 
teo et al. 2002) and free-free emission from the same halos 
that reionize the universe (Oh & Mack 2003), the detec- 
tion and imaging of large-scale structure prior to reioniza- 
tion breakthrough remains a tantalizing possibility within 
range of the next generation of radio arrays. 

The main results of this study can be summarized as 
follows. 

• At epochs when the IGM is still mainly neutral, the 
simulated early galaxy population provides enough 
Lya-c photons to decouple T$ from 7cmb- As in 
the same redshift range the IGM is expected to be 
'warm', the 21 cm line would be seen in emission. 

• The rms temperature fluctuations relative to the 
mean, (ST 2 ) 1 / 2 , increase with decreasing angular 
scale, as variance is larger on smaller scales. The 
signal peaks at an epoch when several high density 
neutral regions are still present, but H n occupies 
roughly half of the volume. 

• Depending on the redshift of reionization break- 
through, broad-beam observations at frequencies < 
150 MHz (below 100 MHz for the 'early reioniza- 
tion' scenario) with the next generation of radio tele- 
scopes should reveal angular fluctuations in the sky 
brightness temperature in the range 5-20 mK (la) 
on scales < 5 arcmin. 
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